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ABSTRACT 

We have studied the O and Ne absorption features in the X-ray spectrum of 
Cyg X-2 observed with the Chandra LETG. The O absorption edge is represented 
by the sum of three absorption-edge components within the hmit of the energy 
resolution and the photon counting statistics. Two of them are due to the atomic 
O; their energies correspond to two distinct spin states of photo-ionized O atoms. 
The remaining edge component is considered to represent compound forms of 
oxide dust grains. Since Cyg X-2 is about 1.4 kpc above the galactic disk, the 
H column densities can be determined by radio (21 cm and CO emission line) 
and observations with relatively small uncertainties. Thus the O abundance 
relative to H can be determined from the absorption edges. We found that 
the dust scattering can affect the apparent depth of the edge of the compound 
forms. We determined the amplitude of the effect, which we consider is the 
largest possible correction factor. The ratio of column densities of O in atomic 
to compound forms and the O total abundance were respectively determined to 
be in the range l.ltH to 2.8+5:5 (ratio), and 0.63 ±0.12 solar to 0.74 ±0.14 solar 
(total), taking into account the uncertainties in the dust-scattering correction 
and in the ionized H column density. We also determined the Ne abundance 
from the absorption edge to be 0.75 ± 0.20 solar. These abundance values are 
smaller than the widely-used solar values but consistent with the latest estimates 
of solar abundance. 
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Subject headings: ISM: abundances — ISM: dust, extinction — X-rays: individ- 
ual(Cyg X-2) — X-rays: ISM 

1. Introduction 

The metal abundance in the interstellar medium (ISM) is an important parameter for 
the understanding of the chemical evolution of the Universe. However, large uncertainties 
remain in our knowledge of this abundance. The solar abundance is sometimes regarded as 
the average of our Galaxy. However, UV absorption lines due to interstellar gas and optical 
absorption lines in the atmosphere of young B-stars have indicated that the metal abundance 
in the ISM of our Galaxy is on average about two thirds of the solar abundance (Savage & 
Sembach 1996; Snow & Witt 1996). However, both methods contain uncertainties because 
the former method is only sensitive to matter in atomic forms, while the abundance of B-star 
atmospheres may be different from the ISM. Indeed, Sofia & Meyer (2001) recently claimed 
that B-star abundances are lower than the ISM abundance and that the solar abundance is 
close to the ISM abundance. Moreover, the O solar abundance itself is not determined well. 
Recent values by Grevesse & Sauval (1998) and Holweger (2001) are significantly lower than 
the values in the widely used table of Anders & Grevesse (1989). 

High resolution X-ray spectroscopy of absorption features due to the interstellar matter 
is a powerful tool to measure the amount of interstellar metal. It is sensitive to both the 
atomic (gas) and compound (molecular and dust grain) forms and can distinguish them by 
the chemical shifts. There are absorption lines and edges of many elements in the X-ray 
energy range. We can determine the column densities of these elements in different forms 
separately, integrating towards the X-ray source. However, the observation requires both 
high energy resolution and good statistics. Schatternburg & Canizares (1986) determined 
the O column density towards the Crab nebula with the transmission grating on the Einstein 
observatory. They could not separate atomic and compound forms. Also they were not able 
to determine the O abundance since the H column density towards the Crab nebula was 
not known. Paerels et al. (2001) analyzed the O and Ne absorption structures of the X- 
ray binary X0614+091 observed with the low energy transmission grating (LETG) on the 
Chandra observatory. They successfully separated the absorption features of O in atomic 
and compound forms. They also determined the abundance ratio of Ne to O. However, they 
could not determine the O abundance because the H column density toward the source is not 
known. Schulz et al. (2002) determined the column densities of O, Ne, Mg, Si, and Fe towards 
Cyg X-1 with the Chandra high energy transmission grating (HETG) data. They discussed 
the abundances of O, Ne, and Fe. However, it was necessary to assume the abundances of 
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all other elements to be solar, in order to estimate the H column density towards Cyg X-1 
from the X-ray spectrum. 

In this paper we analyze the spectrum of Cyg X-2 obtained with the Chandra LETG and 
HRC (High Resolution Camera)-!. Cyg X-2 is a Low-Mass X-ray Binary (LMXB) located 
at {l,b)— (87.33°,— 11.32°) and its distance is determined to be 7.2 kpc from the optical 

counterpart and X-ray burst with mass ejection (Cowley, Crampton, & Hutchings 1979; 
Smale 1998; Orosz & Kuulkers 1999). Smale (1998) claimed a distance of 11.6 kpc, however, 
this is because he assumed the Eddington luminosity for a hydrogen rich atmosphere. If 
we assume a helium rich atmosphere, which is more realistic for type-1 bursts (Lcwin, van 
Paradijs, & Taam 1995), the distance is consistent with 7.2 kpc. Adopting this distance, Cyg 
X-2 is located ~ 1.4 kpc above the Galactic disk. Thus the H column density towards Cyg 
X-2 can be estimated from 21 cm (H I), CO (H2), and Hq, (H II) observations with relatively 
small errors. Because of the featureless X-ray spectrum of the LMXB, the high brightness, 
and the well-determined H column density to the source, Cyg X-2 is one of the best X-ray 
sources for the study of the interstellar abundance with X-ray absorption. 

In the next section, we will first discuss the H column density towards Cyg X-2. Then 
in the following sections we will show the observational results and discuss the O and Ne 
abundances and the uncertainties of the present results. 

Throughout this paper, we quote single parameter errors at the 90 % confidence level 
unless otherwise specified. 

2. H column density towards Cyg X-2 

Interstellar hydrogen is distributed in three different forms; atomic (H I), molecular 
(H2), and ionic (H II). 

The column density of atomic hydrogen is measured by the 21 cm radio emission. In 
the H I map of Dickey & Lockman (1990), the H column density (A^h) is given at grid 
points of 0.8 degree separation. By interpolating between the values of nearby points, the 
column density in the direction of Cyg X-2 is estimated to be A^H(atom)=2.17 x 10^^ cm~^. 
Since the scale height of atomic hydrogen is ~ 140 pc, the total column density of our 
galaxy can be regarded as the column density to Cyg X-2. The standard deviation of the A^h 
values within two degrees of Cyg X-2 is 0.1 x 10^^ cm"^. The power spectra of the small-scale 
spatial fluctuation of the galactic A^H(atom) were measured for different directions (Crovisier 
& Dickey 1983; Dickey et al. 2001). On the spatial scales less than 1 deg, the power-law 
index of the power spectra is in the range from —3 to —4. Assuming the same shape of the 
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power spectrum for the spatial fluctuations in the vicinity of Cyg X-2, we estimate the rms 
amphtude of the spatial fluctuation less than 2 degrees to be 0.12 x 10^^ cm~^ and regard 
this as the 1-a error of the estimation. 

Miramond & Naylor (1995) determined the upper limit of molecular hydrogen column 
density towards Cyg X-2 from the upper limit of CO emission line in the direction of Cyg 
X-2 to be 0.05 x 10^^ cm~^ (90 % significance level). According to them, the contribution 
of H2 in A^H is smaller than 0.10 x 10^^ cm~^. 

The total neutral H column density (A^hi + SAhj) is thus estimated to be 2.17lo;2o ^ 
10^^ cm~^. On the other hand, the reddening towards Cyg X-2 is determined to be E(B — 
V) ~ 0.4 (McChntock et al. 1984). Combining this with the relation, A^hi + 2A^H2 = 5.8 x 
10^^ E{B-V) atoms cm'^ mag"^ (Bohlin, Savage, & Drake 1978), we obtain A'hi + SA'hj ~ 
2.3 X 10^^ cm~^, which is consistent with the above estimation. 

A large fraction (~ 97 %, Mathis 2000) of ionized H atoms exists in the warm ionized 
interstellar medium (WIM) of temperature T ~ 8000 K, distributed in our galaxy with a scale 
height /i of ~ 900 pc and a volume filling factor / of > 0.2. In the WIM, oxygen is effectively 
ionized through the charge exchange process. The number densities of neutral/ionized oxygen 
and hydrogen are related to A^(0 II)/A^ (O I) = (8/9) A^(H II)/A^(H I) (Field & Steigman 
1971). According to more recent model calculations by Sembach et al. (2000), the ionization 
fraction of Ne is also similar to that of H; in their 'COMPOSITE' model, the values of 
log(A^(X*)/A"(Xtotai)) where X — O, or Ne and i represents the ionization states are, 
respectively, -0.73 and -0.09 for H I and H II, -0.71 and -0.09 for O I and O II, and -0.70 
and -0.10 for Ne I and Ne II. Therefore, A^(Xtotai)/A^(Htotai) = N{X^)/N{}i I) is a good 
approximation for gas-phase O and Nc and we do not need to consider ionized H when we 
calculate atomic O and Ne abundances. 

On the other hand, since interstellar dust grains are thought to co-exist with the WIM, 
we need to take into account the column density of ionized H when we estimate the abundance 
of O in compound forms. We estimated it from the emission line intensity in the direction 
of Cyg X-2 obtained from the Wisconsin H-Alpha Mapper (WHAM) survey (Haffner et al. in 
prep) ^. The intensity is 14.6±3.8 Rayleighs. From the line intensity, the total emission mea- 
sure in the column towards Cyg X-2 is estimated to be (1.28±0.33) x 102°(T/8000 K) cm'^ 
where we apphed an extinction correction factor of 1.3 assuming — 3.1 (Mathis 1990). 
Assuming the spatial distribution described above, the emission measure and the column 
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density are respectively related to the local ionized H density in the galactic plane, uq, by, 



and 




where d and b are the distance and the galactic latitude of Cyg X-2, respectively. Prom 
these we obtain A^H(ion) = 0.70 ± 0.18 x 10^^(//0.2)^/^ cm~^. As a nominal value, we 
adopted A^H(ion) = 0.70 ± 0.18 x lO^^ cm-^ (/ = 0.2, Mathis 2000) and thus A^H(total) = 
2.87l°:| X 10^^ cm-2. Even for / = 1, A^H(ion) does not exceed 1.57 ± 0.40 x lO^^ cm-^. 

In Table 1, we summarize the H column densities estimated in this section. 

3. Analysis and results 

3.1. Data reduction 

Cyg X-2 was observed with the LETG/HRC-I for 30 ks on 2000 April 24 (obsID 87). 
We retrieved the archival data from the CXC (Chandra X-ray Center). Cyg X-2 is known 
to show both intensity and spectral variations. However, since the X-ray flux varied only 
~ 20 % during the observation and since we are interested only in the interstellar absorption 
features, we integrated all the data after the standard data screening. 

Since the spectrum in the archival data we retrieved is known to contain some processing 
errors^, we reprocessed the data according to the standard methods described in threads in 
CIAO 2.2.1 (CALDB 2.10). We first obtained the spectra of the positive and negative orders. 
Then we summed them, rebinned to 0.025 A bins, and subtracted the background. The 
background was only ~ 2 % of the source counts in the wavelength region we used in the 
later analysis. 

HRC-I has almost no energy resolution. The spectrum we obtained contains photons not 
only of the first order dispersion but also of the higher orders. We subtracted the higher order 
spectra utilizing the method described in Paerels et al. (2001). The m-th order spectrum of 
photons with a wavelength of A appears in the wavelength bin, A = mX, of the first order 
spectrum. Since the LETG/HRC-I system has no detection efficiency below the wavelength 



^see http://asc.harvard.edu/ciao/threads/spectraJetghrcs/ 
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Ao, where Aq = 1 A, we can neglect the contribution of the m-th order spectrum in the 
wavelength range A < mAo- From the spectrum in the range Aq < A < 2Ao and the efficiency 
ratio of the second to the first order spectra, we can estimate the second order spectrum 
in the range 2Ao < A < 4Ao and subtract it from the total spectrum to obtain the first 
order spectrum in Aq < A < SAq. In the same manner, we can subtract higher order spectra 
from all the wavelength range. We calculated the efficiency ratio from the effective areas of 
the LETG/HRC-S system provided at the CALDB WWW page ^ because the higher order 
efficiencies are not available for the LETG/HRC-I system and because the efficiency ratios 
are not dependent on the detector. The calibration files we used are updated compared to 
those used by Paerels et al. (2001). We show in Figure 1 the spectra with and without higher 
orders. 

We estimated the statistical errors of the first order spectrum obtained with the above 
method considering the error propagations. Because of the error propagations, the statistics 
of the spectral bins are no longer completely independent. However, the increase of the 
statistical errors after the higher order subtraction is smaller than 3% in the wavelength 
ranges we use in the following analysis. We thus treat them to be independent in the 
spectral fittings. 

According to the LETGS calibration report at the CXC ^, the wavelength calibration 
of LETGS spectra is accurate to 0.02 %. This is smaller by an order of magnitude than the 
statistical errors in the wavelengths determined in the following analyses. 

3.2. O absorption structure and the optical depth 

In the observed spectrum of Cyg X-2, we clearly found two absorption lines and an 

absorption edge near A ~ 23 A besides the instrumental O absorption line and edge (see 
Figure 2). The centroid energies of the two absorption lines were determined to be 23.505 ± 
0.007 A and 23.360 ± 0.014 A, and thus they were identified with the ls-2p resonant lines 
of the atomic and compound O, respectively. However, since the absorption lines arc satu- 
rated (r 3> 1), the absorption column densities cannot be determined accurately from the 
absorption lines. We thus focus on the absorption edge. 

The absorption edge structure was extended compared to the instrumental O K-edge, 



3letgs.NOGAP_EA_001031.mod for the 1st order and letgs_EA_001031_om.mod for m-th order, at 
http: / /asc. harvard.edu/cal/Links /Letg/User /Hrc_QE /E A /high_orders / 

"^see http://cxc.harvard.edu/cal/Links/User/calstatus.html 
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which suggests that the edge structure consists of absorption edge components with different 
edge energies. 

A comphcated structure is expected for the absorption edge of an atomic O, firstly be- 
cause photo-ionized O atoms take two distinct final states with different total spin numbers, 
and secondly because there are many absorption lines which cannot be resolved with the 
LETGS near the absorption edges (McLaughlin & Kirby 1998). The theoretical values of 
edge energies calibrated with ground experiments are 544.03 eV (22.790 A, spin 3/2) and 
548.85 cV (22.590 A, spin 1/2) (McLaughlin & Kirby 1998). On the other hand, the edge 
energy of the compound O depends on the chemical state. Thus the absorption edge of the 
compound forms can have an extended structure. The edge energy must be lower than the 
edge energy of atomic O, but higher than the absorption line energy. According to Sevier 
(1979), the edge energy of O in metal oxide is 532.0 ± 0.4 eV (23.31 ± 0.02 A). 

All the O edge structures may not be identified because the statistics and the energy 
resolution of the present data are limited. Thus, we first apply a rather simple model to the 
observed spectrum, then later try more and more complicated models. 

We restricted the wavelength range of model fits in order to avoid the contributions of 
absorption edges of other elements. We did all the analysis described below with several 
different wavelength ranges. We found that the O absorption optical depth determined from 
the analysis is not dependent on the wavelength ranges of the fits, so long as the lower and 
upper hmits of the range are within the ranges 17.5-19.0 A and 27.0-29.0 A, respectively. 
Thus we show here the results with the fitting range of 18.5 A-28.0 A. 

We assumed a power-law function for the continuum spectrum. We also included the ab- 
sorption by neutral matter represented by the tbabs model in the 'sherpa' program (Wilms, 
Allen, & McCray 2000) with the O abundance fixed to 0. We included an O absorption 
edge with a separate absorption edge component represented by edge model in the sherpa 
program. We fixed the H column density of the tbabs absorption model to 2.17 x 10^^ cm~^ 
and fixed the helium abundance to the solar value. We performed the fits with the abun- 
dances of other atoms fixed at several different values between 0.5 and 1.5 solar. We found 
that the dependence of the best-fit O absorption optical depths on the abundance of other 
elements is only about ten percent of the statistical errors. Thus we only show the results 
with the abundance of other elements fixed at the solar values. 

In the first model, we included two oxygen edge components. Wc treated the edge ener- 
gies and the optical depths of the two 'edge' components as free parameters. We represented 
the two O absorption lines with negative Gaussian functions. We fixed the centroid wave- 
lengths to 23.505 A and 23.360 A, and the widths to 0.04 A. The free parameters of the fits 
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were the two O edge energies and optical depths, the normahzation factor and the power-law 
index of the continuum, and the normalization factors of the two Gaussian absorption lines. 

We show the results of the fits in the first data column of Table 2. The best fit energies 
of the two edges were 22.76ig:J^ ^ and 23.1llg:g| A, respectively. The former is consistent 
with the atomic O edge of the spin 3/2 final state within the statistical error. On the other 
hand, the latter is in the middle of the atomic and the instrumental O absorption edges. 
Thus, it is likely to represent an O edge of compound forms. 

Next, in order to investigate possible complex structures of the edges and the model 
dependence of the absorption optical depths, we increased the number of edge components 
one by one. In the fits, we allowed both the energies and optical depths of all the edges to 
vary. In the second and third data columns of Table 2, the results are summarized. The 
value shows a significant improvement from the double to the three edge models (99.6 % 
confidence with an F-test). The edge energy of the third edge is 22.58 A= 549.1 eV. This 
energy is close to the edge energy for a final spin 1/2. In the raw spectrum (Figure 2), we 
can find this edge structure. On the other hand, the value did not improve from the 
three-edge to the four-edge models. Thus we conclude that the observed O edge structure 
can be represented with the three-edge model within the limit of the present statistics and 
the energy resolution; two of the edges represent absorption by atomic O corresponding to 
the two different O-ion spin states, and the remaining one represents absorption by O in 
compound forms. 

Finally, in order to determine the spectral parameters and their statistical errors for the 

three-edge model more precisely, we performed the spectral fit with the three-edge model 
with the two edge energies of the atomic O fixed at the theoretical values. In the last column 
of Table 2, we show the results. We estimated the statistical error for the total optical depth 
for the atomic O edge from the contour map to find the total optical depth of 0.43 ±0.14. 
The optical depth of compound forms is 0.26 ±0.11. 

To estimate the column density of O from the absorption optical depth, we adopted 
the absorption cross section of Verner & Yakovlev (1995): dabs = 4.98 x 10~^^ cm^ both for 
atomic and compound forms. We obtained iVo,atom = 8.6 ± 2.8 x 10^'' cm~^ and A^o.comp = 
5.2 ±2.2 X 10^^ cm-2. 

Assuming the A^h values of 2.17lo:2o x 10^^ cm"^ for atomic O and l.m'^^^f^ x 10^^ cm"^ 
for compound forms, we estimate the abundance of O to be 0.47 ± 0.16 solar (atomic), 
0.21 ± 0.09 solar (compound), and 0.68 ± 0.13 solar (total), where we assumed the O solar 
abundance of 8.51 x 10~^ (Anders & Grevesse 1989). The error of the total abundance was 
estimated from the contour map of the spectral fit, taking into account the difference of 
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A^H values for gas and compound forms and their systematic errors. 

3.3. Absorption structure of Ne 

In the observed spectrum of Cyg X-2, we also clearly found an absorption edge at 14.3 
A, which is identified with the Ne K-edge (Figure 3). We performed spectral fits to the 
wavelength range of 12.0-15.5 A. In addition to the edge, there is an absorption line at 14.6 
A. Because the hne wavelength is close to the edge, we included an absorption hne in the 
model function. We employed the model function similar to that used for the O features. 
We represented the Ne edge with a single edge model. 

In Table 3, we show the results of the fit. The absorption edge wavelength was consistent 
with that of the neutral Ne (14.25 A; Verner & Yakovlev 1995). Prom the optical depth 
determined from the fit and the theoretical absorption cross section of Ne I (Verner & 
Yakovlev 1995), we determined the column density and abundance of Ne to be 2.00 ± 0.51 x 
10^^ cm~^ and 0.75 ± 0.20 times the solar abundance of Anders & Grevesse (1989). 

Although the absorption-line energy, 14.599±0.012 A, is close to the resonant line energy 
of Ne II (14.631 A, Behar & Netzer 2002), it is outside the 90% statistical error domain. 
Thus we can not unequivocally identify this absorption line. However, if we assume it is 
Ne II, the column density estimated from the line is 8.0 ± 2.4 x 10^^ cm~^, which is about 
1/3 of Ne I. Thus if the absorption line is due to Ne II, Ne II can be attributed to the WIM. 

4. Discussion 

We have analyzed the O and Ne absorption structures in the X-ray spectrum of Cyg X-2 
observed with the Chandra LETG/HRC-1. The O edge is represented by the sum of three 
edge components within the limit of present statistics and energy resolution. Two edge 
components represent absorptions by atomic O; their edge energies are consistent with the 
theoretical values corresponding to the two distinct final O states of different spin numbers. 
The edge energy of the remaining edge is lower than the other two, and is likely to represent 
O in compound forms. Prom the depth of the edges, we estimated the O column densities 
with atomic and compound forms separately. 

We tried a four-edge component model in the spectral fits. We find the sum of the 
optical depths of all the edge components does not change significantly from three to four 
edge models. Thus we conclude that the total O column density does not depend significantly 
on the number of edge components assumed in the spectral models. We estimated the atomic. 
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compound, and total O abundances to be 0.47 ± 0.16, 0.21 ± 0.09, and 0.68 ± 0.13 times 

the solar abundance of Anders & Grevcssc (1989) respectively. The Ne edge is represented 
with a single edge model and the Ne abundance is estimated to be 0.75 ± 0.20 solar. In the 
errors of abundance quoted above, we included the statistical errors of the spectral fits, and 
the systematic errors of A^h estimated from the spatial fluctuation of 21cm and emissions 
and the upper limit of CO emission. 

There is an additional systematic error in the abundance of O in compound forms, arising 
from the uncertainty in the spatial distribution of the WIM. If we consider the extreme case 
in which its volume filling factor is 1, the O abundance reduces to 0.16 ± 0.07 (compound 
form) and 0.63 ± 0.12 (total). 

It is very unlikely that the O and Ne absorption edges are associated with the circum- 
stellar matter in the Cygnus X-2 binary system. Suppose that the absorption column of 
iVn ~ 2 X lO^^cm"^ is located in a distance r from the X-ray source. The ionization param- 
eter, ^ = Lx/ (nr^) ~ Lx/ {Nar) must be smaller than ~ 10 so that O is neutral (Kallman & 
McCray 1982). With Lx — 10^^ erg sec~^, we find r > 10^^ cm. This is much larger than the 
size of the binary system (~ 10^^ cm, Cowley, Crampton, & Hutchings 1979) This is consis- 
tent with the facts that all emission and absorption line features found in LMXBs are from 
highly ionized atoms (e.g., Asai ct al. 2000; Cottam et al. 2001) and that low-ionization emis- 
sion lines from the Cyg X-2 system are from near the companion star (Cowley, Crampton, 
& Hutchings 1979). 

About one third of O contributing to the absorption edges is in compound form. Dust 
grains are considered to contribute to this component. In such a situation, the optical depth 
for dust scattering of X-rays is not negligible compared to the absorption if we assume 
a typical dust radius of ~ 0.1 fira. Scattered X-rays form an extended halo on the few 
arcminute scale (Hayakawa 1970; Mauche & Gorenstein 1986; Prcdehl & Schmitt 1995). 
Since the scattering cross section is dependent on the X-ray energy, the energy spectrum of 
the unscattered X-rays is modified by scattering. However, if multiple scattering is negligible 
and if the spatial distribution of dust is uniform over the spatial scale corresponding to the 
dust halo, the energy spectrum of total photons, i.e., the sum of the unscattered core and 
the scattered halo, is not modified by dust scattering. This is not the case for the present 
energy spectrum obtained with the LETG, since we derived the energy spectrum only from 
the central ~ 1 arc second of the X-ray image. Because the scattering cross sections show 
anomalous features around absorption edges (Mitsuda et al. 1990), the dust scattering will 
affect the estimation of the absorption column density; the scattering cross section is smaller 
at the energy just above the edge (at the wavelength shorter than the edge), while the 
absorption cross section is larger. This leads to an underestimation of the O absorption 



column density. We estimated the amplitude of this effect in the following way. 

The total scattering cross section of an X-ray photon of a wavelength A with a dust 
grain of a radius a is given by 



under the Rayleigh-Gans approximation (van de Hulst 1957), where N is the number of 

atoms per unit volume, Tq is the classical electron radius, and /i and /2 are the atomic 
scattering factors (Henke, Gullikson, & Davis 1993). The Rayleigh-Gans approximation is 



Because the factor |/i-|-i/2| becomes small near the absorption edge, the dust scattering 
cross section is reduced at the edge. Assuming typical chemical compositions and densities 
of dust grains (FeSi04, FeSiOs, Mg2Si04 and MgSiOa), we estimate the reduction of the 
scattering cross section at the O edge to be in the range. 



Assuming the number density of dust to be proportional to oc a~^'^ with a cut off at Omax 
(Mathis, Rumpl, & Nordsieck 1977), we calculate the reduction of the scattering cross section 
per O atom. 



Thus the reduction of scattering cross section can be a few tens of percent of the absorption 
cross section, dabs = 4.98 x 10~^^ cm^. The maximum grain radius contributing to X-ray 
scatterings has been estimated from the spatial size of dust scattering halos. Mauche & 
Gorenstein (1986) and Mitsuda et al. (1990) suggested Omax = 0.05 — 0.1 /xm, but recently 
Witt, Smith, & Dwek (2001) suggested the presence of grains with radii as large as 0.5/xm. 
Adopting Omax = 0.1 //m, and assuming that O of compound forms are all in dust grains, 
the column density of O of compound forms increases from 5.2 x 10^^ cm~^ to (6.1 — 7.3) x 
10^^ cm~^. At the O edge energy, the Rayleigh-Gans approximation is valid only for Cmax 
0.1 jim. and the scattering cross section saturates at flmax ^ 0.1 yum (Alcock & Hatchett 1978; 
Smith & Dwek 1998). Thus the correction factor does not increase for larger grains and the 
correction factor we applied above can be regarded as the maximum correction. 

The O abundance values calculated for four different cases (with and without correction 
for dust scattering, and two different values of the volume filling factor of the WIM) are 
summarized in Table 4. The total O abundance is between 0.63 ±0.12 and 0.74 ±0.14 times 




valid for (27ra/A)|iVroA2 



/(27r)(/i + i/2)|«l. 
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solar. On the other hand the Ne abundance is 0.75 ± 0.20. Our results are more consistent 
with the recent solar abundances than the most widely used 'old' values. For example , the 
O and Ne solar abundances by Holweger (2001) are 0.65 and 0.81 times the values in Anders 
& Grevesse (1989), respectively. 

Another important parameter we can estimate from the present X-ray observation is 
the gas to dust ratio of the O column densities, which is estimated to be 2.2^^'^, 2.8^^5, 
1-71q 9, and 2.2ti'? for cases 1 to 4 of Table 4, respectively. This should be compared with 
the value estimated from the gas abundance from interstellar UV absorption lines and the 
assumption that the total abundance is solar. Adopting the O gas abundance of Cartlcdge, 
Meyer, & Lauroesch (2001) and solar abundances by Anders & Grevesse (1989) and Holweger 
(2001), we obtain 0.7 and 1.7, respectively. Our result is again consistent with the latest 
solar abundance value. 

For O, the present X-ray result contains an additional uncertainty due to the correc- 
tion for dust scattering. We can avoid this if we include the spectrum from the scattering 
halo in the analysis. This requires non-dispersive high resolution spectrometers such as 
microcalorimeters, which will be used on board future X-ray missions, such as Astro- E2. 

We thank J. Ichimura and K. Masai for discussions on the O and Ne absorption structure, 
P. Edwards for his careful review of the manuscript. We appreciate the Hq, survey of the 
Wisconsin H- Alpha Mapper funded by the National Science Foundation. We are also grateful 
to the anonymous referee for his valuable comments, in particular, suggesting the importance 
of H II in the H column density. 
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Table 1. Hydrogen column density to Cyg X-2 



Form of H H column density 



HI 2.17 ± 0.20 X 10^1 cm-2 

H2 < 0.10 X 10^1 cm-2 

H II (nominal) 0.70 ± 0.18 x lO^i cm'^ 

H II (maximum )(^) 1.57 ± 0.40 X 10^1 cm-2 

H I+H2 2.17l°|2 X 1021 cm-2 

All forms (nominal) 2.87i[J:| x 10^^ cm-^ 

All forms (maximum) 3.74 ± 0.45 x lO^^ cm-^ 



Column density to be used for the calculation of O in the atomic form and Ne. 
*^2^ Column density to be used for the calculation of O in the compound forms. 
(^^Volume filling factor of 0.2 is assumed for WIM. 
(^^Volume filling factor of 1 is assumed for WIM as an extreme case. 



Table 2. Results of spectral fits of O K-edge absorption features 



Parameters 



Two Edges Three Edges Four Edges'^^ Three with two fixed^^^ 



Power-law 
normalization 
photon index 



(photons s-i keV-i cm-2@555.2 eV) 2.52 ± 0.07 2.53 ± 0.07 



0.74 ±0.17 



0.69 



+0.18 



2.53 
0.70 



2.53 ± 0.07 
0.71 ±0.17 



Absorption line 1 
Centroid 

Equivalent width 

Absorption line 2 
Centroid 

Equivalent width 

Absorption Edge 1 
Edge energy 

Optical depth 

Absorption Edge 2 
Edge energy 

Optical depth 

Absorption Edge 3 
Edge energy 

Optical depth 

Absorption Edge 4 
Edge energy 

Optical depth 

X^/dof 



(eV) 

(A) 

(eV) 



(eV) 

(A) 

(eV) 



(eV) 

(A) 



(eV) 

(A) 



(eV) 
(A) 



(eV) 
(A) 



527.48 
(23.51) 
1.41 ± 0.24 



530.76 
(23.36) 
0.97 ±0.28 



536.4t°| 
(23.1ltH«) 
0.27 ±0.10 



(22, 



7f;+0.i3-, 



41+°10 



527.48 
(23.51) 
1.41 ±0.24 



530.76 
(23.36) 
0.97 ±0.28 



536.2(3) 
(23.12) 
0.24 ±0.12 



542.3(3) 

(22.86) 
(, 24+0-29 



549.1(3) 

(22.58) 
22+°-^3 



527.48 
(23.51) 
1.41 



530.76 
(23.36) 
1.0 



536.2 

(23.13) 
0.24 



542.3 
(22.86) 
0.19 



544.9 
(22.75) 
0.09 



549.2 
(22.58) 
0.18 



527.48 
(23.51) 
1.41 ±0.24 



530.76 
(23.36) 
0.97 ±0.28 



536.1 ±1.5 

(23.13 ±0.06) 
0.26 ±0.11 



544.03 
(22.79) 
0.23; 



+0.16 
0.11 



548.85 
(22.59) 

o.2ot°;lf 



400.45/371 397.40/369 397.163/367 



398.21/371 



(^) Statistical errors are not estimated for this model. 

(■^)Edgc energies of atomic O are fixed at the theoretical values. 

(3) Statistical errors are not estimated because the edge energies are strongly coupled. 
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Wavelength (A) 

Fig. 1. — The Cyg X-2 spectrum obtained with the Chandra LETG/HRC-I. The raw LETG 
spectrum containing higher order photons is shown with a thin hne, while the spectrum after 
higher-order subtraction (see text) is shown with a thick hne. 



Table 3. Results of spectral fits of Ne K-Edge 



Parameters 




Best fit value 


Power-law 






normalization 


(photons s~i keV~i cm~2@lkeV) 


1.685 ± 0.023 


photon index 




1 QQQ + O.OTO 
1.0^0_Q gy2 


Absorption line 






Centroid 


(eV) 


849.24 ± 0.68 




(A) 


(14.559 ±0.012) 


Equivalent width 


(eV) 


0.541 ±0.160 


Absorption Edge 






Edge energy 


(eV) 


866 ±5 




(A) 


(14.311HD 


Optical depth 




0.059 ± 0.015 






149.724/133 
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Fig. 2. — Spectral fit in the wavelength range including the O absorption structures. The 
whole wavelength range used in the spectral fits is shown in the lower graph, while in the 
upper graph the narrow wavelength range containing O edges and absorption lines is ex- 
panded. In the upper graph and in the upper panel of the lower graph, the observed 1st 
order spectrum is shown with data points with vertical bars which represent \-a statistical 
errors. The best-fit model function of the three edge-component model (the last column of 
Table 2) convolved with the telescope and the detector response functions is shown with a 
thick line. The model function without O edges and absorption lines is shown with a thin 
line, in order to show the instrumental O absorption structure. Two absorption lines and an 
extended O edge clearly exist. In the lower panel of the lower graph, the residuals of the fit 
are shown 
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Wavelength (A) 



Fig. 3. — Spectral fit in the wavelength range including the Ne absorption structures. The 
observed 1st order spectrum with 1-a error bars and the best fit model function convolved 
with telescope and detector response functions are shown in the upper panel, and the resid- 
uals of the fit are shown in the lower panel. 
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Table 4. Summary of interstellar O and Ne abundances 



Atom 


Form 


Case l'^ 


Case 2^ 


Case 3*^ 


Case 4^1 





gas phase 


0.47 ±0.16 










compound 


0.21 ± 0.09 


0.16 ±0.07 


0.27 ±0.12 


0.21 ± 0.09 




total 


0.66 ±0.12 


0.63 ±0.12 


0.74 ±0.14 


0.68 ±0.14 


Ne 




0.75 ± 0.20 









Note. — All abundance values are in units of the solar values of Anders & Grevesse (1989). 

^Case 1: Volume fiUing factor of 0.2 is assumed for the WIM. No correction is made for 
dust scattering. 

''Case 2: Volume filling factor of 1.0 is assumed for the WIM as an extreme case. No 
correction is made for dust scattering. 

'^Case 3: Volume filling factor of 0.2 is assumed for the WIM. Correction which is considered 
to be the maximum is made for dust scattering. 

'^Case 4: Volume filling factor of 1.0 is assumed for the WIM as the extreme case. Correc- 
tion which is considered to be the maximum is made for dust scattering. 



